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Abstract 

Background: How plants respond to dehydration stress has been extensively researched. However, how plants respond 
to multiple consecutive stresses is virtually unknown. Pre-exposure to various abiotic stresses (including dehydration) may 
alter plants' subsequent responses by improving resistance to future exposures. These observations have led to the 
concept of 'stress memory' implying that during subsequent exposures plants provide responses that are different from 
those during their first encounter with the stress. Genes that provide altered responses in a subsequent stress define the 
'memory genes' category; genes responding similarly to each stress form the 'non-memory' category. 

Results: Using a genome-wide RNA-Seq approach we determine the transcriptional responses of Arabidopsis plants that 
have experienced multiple exposures to dehydration stress and compare them with the transcriptional behavior of 
plants encountering the stress for the first time. The major contribution of this study is the revealed existence of 
four distinct, previously unknown, transcription memory response patterns of dehydration stress genes in A.thaliana. 
The biological relevance for each of the four memory types is considered in the context of four overlapping strategies 
employed by a plant to improve its stress tolerance and/or survival: 1) increased synthesis of protective, 
damage-repairing, and detoxifying functions; 2) coordinating photosynthesis and growth under repetitive stress; 
3) re-adjusting osmotic and ionic equilibrium to maintain homeostasis; and 4) re-adjusting interactions between 
dehydration and other stress/hormone regulated pathways. 

Conclusions: The results reveal the unknown, hitherto, existence of four distinct transcription memory 
response types in a plant and provide genome-wide characterization of memory and non-memory dehydration 
stress response genes in A.thaliana. The transcriptional responses during repeated exposures to stress are 
different from known responses occurring during a single exposure. GO analyses of encoded proteins suggested 
implications for the cellular/organismal protective, adaptive, and survival functions encoded by the memory genes. The 
results add a new dimension to our understanding of plants' responses to dehydration stress and to current models for 
interactions between different signaling systems when adjusting to repeated spells of water deficits. 



Background 

Environmental stresses exert evolutionary pressure on 
organisms, which, in turn, have developed sophisticated 
responses to cope and to survive. Plants react to unfavor- 
able conditions by dynamically changing physiological 
behavior and expression levels of implicated genes. Pre- 
exposure to various abiotic stresses (salinity, cold, high 
temperature) may alter plants' subsequent responses by im- 
proving resistance to future exposures [1,2]. Pre-treatment 
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(priming) with hormones (jasmonic acid, J A, abscisic acid, 
ABA, salicylic acid, SA) increased the systemic immunity 
and induced stronger responses from responding genes 
upon subsequent treatments, relative to non-primed plants 
[3-7]. These observations have led to the concept of stress 
memory; implying that during subsequent exposures plants 
provide responses that are different from their responses 
during the first encounter with the stress [8]. 

Drought-triggered dehydration stress is one of the most 
common environmental stresses endured by plants and 
their responses to dehydration stress are extensively re- 
searched at organismal, cellular, and genome levels. Re- 
cently, we demonstrated that Arabidopsis plants subjected 
to several cycles of dehydration/water recovery treatments 
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maintained higher relative water content than plants ex- 
periencing dehydration stress for the first time [9]. Dehy- 
dration stress memory' also affected gene expression from 
response genes. Analysis of 14 dehydration response genes 
after multiple stress treatments revealed the existence of 
two different types of transcriptional responses: genes that 
produced transcripts at similar levels during each stress 
and genes that in a subsequent response significantly in- 
creased their transcript levels and transcription rates rela- 
tive to the first response [9]. Accordingly, genes that 
provide altered responses in a subsequent stress were 
referred to as 'memory genes' to distinguish them from 
'non-memory' genes responding similarly to each stress. 

Here, we explore how broadly memory genes are dis- 
tributed within the dehydration stress responsive fraction 
of the Arabidopsis thaliana genome. Comprehensive 
RNA-Seq analyses of the transcriptomes of plants prior 
to stress, during a first stress, and during a third dehy- 
dration stress exposure revealed an unexpected diversity 
of memory-type responses at the transcriptional level. 
Analyses of the memory genes by GO functional cat- 
egories indicated a biased functional distribution among 
the memory response patterns. Results are interpreted 
in the context of the possible biological relevance of 
these memory responses. 

Methods 

Plant growth and treatments 

Arabidopsis thaliana (Col-0) plants were grown in potting 
soil in growth rooms at 22°C with a 12-h light photo- 
period and light intensity of 180 (imol m" 2 s~\ Repeated 
dehydration stresses were performed by air-drying for 2 h 
followed by a 22 h period of full re-hydration recovery as 
described [9]. RNA-Seq analyses were performed on ros- 
ette leaves from pre-stressed (W) plants, from plants ex- 
posed to the first dehydration treatment (SI), and during 
a third stress (S3) following two stress/recovery cycles. 
Plants from two independent biological samples were 
used. After removal from soil and before initiating the 
stress/recovery cycle, experimental plants were condi- 
tioned overnight in humid chambers and used to establish 
the basal transcript levels (W) for the transcriptome ana- 
lyses. The transcriptional behavior of the marker memory 
genes RD29B and RAB18 were monitored in leaf samples 
taken from fresh plants in soil, from overnight condi- 
tioned plants after their removal from soil, and from 
plants kept under watered conditions throughout the 
stress/recovery as indicators of the normal dehydration 
stress behavior of plants in SI and S2 and internal con- 
trols for the whole genome transcriptome analyses. 

RNA extraction and RNA-Seq library construction 

Leaf tissues were collected and immediately frozen in 
liquid nitrogen. Total RNA was extracted with Trizol 



(Invitrogen Inc. Carlsbad, CA, USA), treated with DNase 
I (Qiagen, Valencia, CA), and purified using Qiagen 
RNeasy Mini Kit. RNA integrity was confirmed on a 
Bioanalyzer 2100 using Nano 6000 LabChip (Agilent 
Technologies, Santa Clara, CA). Complementary DNA 
sequencing library was prepared from the total RNA 
using the mRNA-Seq Sample Preparation Kit (Illumina, 
San Diego, CA). Briefly, poly-adenylated RNA was iso- 
lated from 10 \xg total RNA by Sera-Mag Magnetic 
Oligo-dT beads (ThermoFisher Scientific, Waltham, 
MA). RNA Purified mRNA was fragmented, annealed to 
high concentrations of random hexamers, and reverse 
transcribed. Following second strand cDNA synthesis, 
end repair, and A-tailing, Oligo adapters complementary 
to sequencing primers were ligated to cDNA fragment 
ends. Resultant cDNA libraries were size fractionated on 
an agarose gel, 200 bp fragments excised, and amplified 
by 15 cycles of polymerase chain reaction. Clusters were 
generated from the cDNA sequencing library on the sur- 
face of a flowcell in the Cluster Station (Illumina) by so- 
called bridge amplification. Replicates for the watered, 
SI and S3 sample libraries were each run on a single 
lane in a flowcell on an Illumina GAIIx at the Genomics 
Core Facility at the University of Nebraska-Lincoln. 

Reverse transcription and real-time PCR 

Total RNA isolation and reverse transcription with oligo 
(dT) (18418-012, Invitrogen) were performed as de- 
scribed previously [10]. The amounts of individual genes 
were measured with gene-specific primers by real-time 
PCR analysis with a CyclerlQ real-time PCR Instrument 
(Bio-Rad) and SYBR Green mixture (Bio-Rad). The rela- 
tive expression or amount of specific genes was quanti- 
tated with the 2" AACt calculation [11], according to the 
manufacturer's software (Bio-Rad), where the reference 
gene was ubiquitin. Primers used in real-time RT-PCR 
are in Additional file 1: Table SI. 

Bioinformatics analysis 

Transcriptome sequencing of the watered, SI, and S3 
samples yielded a total of 53.4, 76.8, and 78.4 million 
reads, respectively, summed over the two biological rep- 
licates per sample (see Additional file 2: Table S2). The 
read length for SI and S3 is 75 bases, while for watered 
it is 101 bases. To determine the quality of the replicates 
we performed a least-square simple linear regression for 
each of the three samples. We calculated the R 2 statistic 
(0.96 < R 2 < 0.99) and slope (1.07 < b < 1.09), which provide 
measures of goodness-of-fit and correlation, respectively, 
using the regress function in MATLAB® (version 8.0.0.783 
[R2012b]; The MathWorks™) (Additional file 3: Table S3). 
For use in all further analyses, the Arabidopsis thaliana 
genome and gene models were downloaded from Alumina's 
iGenomes (Ensembl, TAIR10), and gene functional 
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descriptions and Gene Ontology assignments from the 
Arabidopsis Information Resource (release TAIR10, [12]). 

The bowtie (version 2.1.0; [13]) and tophat (version 
2.0.8; [14]) packages were used with default parameters 
to map the RNA sequence reads from watered, SI, and 
S3 to the genome (see Additional file 3: Table S3), and 
to determine the expression quantity of known tran- 
scripts in each sample. The cuffdiff tool from the cuf- 
flinks package (version 2.0.2; [15]) was run with default 
parameters to calculate expression changes and associ- 
ated q-values (False Discovery Rate adjusted p-values) 
for each gene, between the samples SI and water, and S3 
and SI. We further classify genes as being significantly 
differentially expressed when all three of the following 
conditions are met: q < 0.05; | log 2 (fold change) | > 1; 
and the FPKM-normalized expression value of at least 
one sample out of the two needs to be larger than the 
25 th percentile. The output files of cuffdiff 'are further an- 
notated (in-house Perl script) by adding gene functional 
descriptions and GO classifications, and merged into a 
master file containing all data for SI versus water and 
S3 versus SI (Additional file 3: Table S3, Additional file 
4: Table S4, and Additional file 5: Table S5). 

From that master file we determined the 6579 signifi- 
cant drought-responsive genes (SI versus water), and 
using that initial set we then looked at the significant re- 
sponses in S3 versus SI (1963 genes) and the remainder 
that did not respond (4616 genes). We assigned simple 
classifications to the types of response during the first 
stress (+ or -) and the second stress (+, -, or =), combin- 
ing them into six classes: [++], [ — ], [+-], [-+], [+=], 
and [-=] (Table 1; Table 2; Additional file 4: Table S4, 
Additional file 5: Table S5). Two additional classes [=/+] 



Table 1 Dehydrations response and memory genes in 



Arabidopsis thaliana 





Total 


Up 


Down 


Total genes 


33,555 






Dehydration response 


6,579 


3396 


3183 


Memory genes 


1963 


1219 


744 


[+/+] W<S1 <S3 


362 






[-/-] W>S1 >S3 


310 






[+/-] W<S1 >S3 


857 






[-/+] W>S1 <S3 


434 






Non-memory genes 








[+/=] W<S1 =S3 


2,177 






[-/=] W>S1 =S3 


2,439 






Late-response genes 








[=/+] W = SKS3 


798 






[=/-] W = S1 >S3 


573 







and [=/-] contain genes that in SI did not change signifi- 
cantly expression (according to the three criteria for sig- 
nificance, above) compared to pre-stressed levels in W but 
significantly changed transcription in S2. Formally, these 
genes do not belong to the SI dehydration-stress respond- 
ing fraction. 

The raw transcriptome sequence files for watered, SI, 
and S3 have been uploaded, together with gene expression 
result files, to NCBIs Gene Expression Omnibus under 
sequence number GSE48235. 

Results 

Dehydration stress memory response genes of 
Arabidopsis thaliana 

Watered non-stressed plants and plants subjected to one 
or three dehydration stresses by exposure to dry air for 2 
hours followed by 22 h periods of watered recovery in- 
tervals [9] were analyzed. The levels of mRNAs from the 
leaves of these plants were determined by RNA-Seq. 
Transcripts from 33,555 A. thaliana genes were identi- 
fied (Table 1; Additional file 3: Table S3). Of these, 6579 
genes (-20% of the genes in the genome) significantly 
increased or decreased their transcript levels during the 
first dehydration stress (SI), compared to their transcript 
levels in watered plants (W), representing the dehydra- 
tion stress-responsive genes of A. thaliana. Of the genes 
responding in SI, 3396 genes were up-regulated and 
3183 were down-regulated (Table 1; Additional file 4: 
Table S4). A comparison of the transcript levels during 
the first (SI) and third dehydration stress (S3) identified 
1,963 dehydration responsive genes displaying signifi- 
cantly different amounts of transcripts in S3 versus SI. 
Our operational criterion for transcriptional memory is 
that the transcriptional responses to similar stress condi- 
tions must be different. Accordingly, these 1,963 genes 
display dehydration stress transcriptional memory (Table 1; 
Additional file 4: Table S4). 

Four distinct transcriptional memory response-patterns 
were recognized within the SI -responding fraction. We 
designate these transcription patterns as [+/+], [-/-], 
[+/-] and [-/+] memory responses, where the first sign 
indicates higher (+) or lower (-) transcript levels in SI 
relative to the levels in pre-stressed watered plants (W). 
The second sign indicates higher (+) or lower (-) tran- 
scripts in S3 relative to the levels in SI. Non-memory 
genes are denoted as [+/=] or [-/=] indicating that up- 
regulated or down-regulated transcript levels occurring in 
SI, respectively, have similar levels in S3. In addition, 798 
genes were up-regulated in S2 ([=/+] W = S1<S3) and 
573 were down- regulated ([=/-] W=S1>S3). These de- 
layed response patterns are different from the memory re- 
sponses represented by the four SI -responding memory 
types. Here, we focus mainly on the transcriptional pat- 
terns of the Sl-responding gene fraction of Arabidopsis. 
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Table 2 Distribution of dehydration stress responding genes according to GO function 





[+/+] 


[-/-] 


r+/-i 

L ' / J 


[-/+] 


[+/=] 


[-/=] 


[=/+] 


[=/-] 




362 


310 


857 


434 


2177 


2439 


798 


573 


Ribosomal and protein synthesis 


1 


31(10%) 


1 


0 


10 (0.5%) 


135 (5%) 


2 


26 (5%) 


Response to salt 


25 (7%) 


8 (3%) 


84 (10%) 


14 (3%) 


1 54 (7%) 


55 (2%) 


35 (4%) 


24 (4%) 


Response to cold/heat 


28 (8%) 


1 8 (8%) 


53 (6%) 


1 9 (4%) 


147 (7%) 


72 (3%) 


40 (5%) 


46 (7%) 


Response to ABA 


24 (7%) 


6 (2%) 


80 (9%) 


12 (1%) 


143 (7%) 


38 (1%) 


39 (5%) 


24 (4%) 


LEA 


12 (3%) 


0 


4 


2 


9 


2 


5 (0.7%) 


0 


Response to light 


18 (5%) 


25 (8%) 


28 (3%) 


24(6%) 


68 (3%) 


101 (4%) 


47 (6%) 


27 (5%) 


UV 


3 


2 


4 


5 


10 


18 


12 


6 


Intensity 


11 


10 


13 


11 


40 


66 


27 


14 


Red/far red 


2 


12 


7 


3 


14 


16 


6 


7 


Blue 


2 


1 


3 


5 


2 


10 


2 


0 


Orcadian rhythm 


1 


1 


13 


6 


11 


23 


8 


12 


Chloroplast 


12 (3%) 


75 (24%) 


20 (2%) 


22 (5%) 


48 (2%) 


203 (8%) 


26 (3%) 


75 (15%) 


Thylakoid membrane 


0 


53 (17%) 


6 (0.7%) 


9 (2%) 


14 


1 77 (7%) 


7 (1%) 


22 (4%) 


Membrane 


30 (8%) 


23 (7%) 


106 (12%) 


41 (10%) 


312 (14%) 


233 (10%) 


73 (9%) 


71 (14%) 


TM transport, porins 


20 (5%) 


6 (2%) 


62 (7%) 


26 (5%) 


138 (7%) 


95 (4%) 


50 (6%) 


24 (4%) 


Wall/PM 


13 (3%) 


25 (8%) 


30 (3%) 


22 (5%) 


48 (2%) 


129 (5%) 


38 (4%) 


37 (6%) 


Kinases, recepors, signal trans-duction 


4 


16 (5%) 


52 (6%) 


31 (7%) 


1 07 (5%) 


111 (5%) 


25 (3%) 


21 (4%) 


Response to auxin 


9 (2%) 


5 (2%) 


43 (4%) 


10 (1%) 


66 (3%) 


48 (2%) 


12 (1%) 


1 2 (2%) 


Response to ethylene 


4 (1%) 


3 (1%) 


53 (6%) 


8 (2%) 


87 (4%) 


24 (1%) 


12 (%) 


8 (1%) 


Response to GA 


2 (1%) 


5 (2%) 


10 (1%) 


7 (1%) 


13 


18 (0.7%) 


13 (1%) 


5 (1%) 


Response to JA 


7 (2%) 


16 (5%) 


121 (14%) 


1 9 (2%) 


89 (4%) 


46 (2%) 


13 (1%) 


21 (4%) 


Response to SA 


4 (1%) 


5 (2%) 


41 (5%) 


13 (3%) 


63 (3%) 


24 (1%) 


14 (1%) 


12 (2%) 


Transcription factors 


29 (8%) 


6 (2%) 


73 (7%) 


25 (6%) 


161 (7%) 


1 00 (4%) 


51(2%) 


22 (4%) 


AP2/ERF 


5 


1 


16 


2 


19 


10 


1 


3 


bHLH 


2 


1 


16 


2 


11 


16 


5 


0 


homeo_ZIP 


3 


1 


6 


3 


25 


16 


8 


3 


MYB 


4 


1 


8 


3 


16 


9 


10 


4 


ZF 


3 


2 


6 


11 


41 


28 


12 


4 


B_ZIP 


4 


0 


4 


2 


13 


1 


3 


0 


NAC 


1 


0 


7 


1 


13 


1 


3 


0 


GRAS 


1 


0 


2 


0 


3 


6 


1 


2 


HSF 


1 


0 


3 


0 


4 


0 


2 


1 


CCAAT 


3 


0 


0 


0 


2 


3 


2 


0 


WRKY 


2 


0 


5 


1 


14 


4 


1 


1 


Number of genes and percentages per memory group are reported. Only percentages higher than 0.7 are reported; higher than 1% are rounded to the 



nearest integer. 



Of the four different memory classes, 362 genes de- 
noted as [+/+] memory genes were induced in SI and 
induced to higher levels in S3 (Table 1; Figure 1A). All 7 
memory genes described in our earlier study [9] were 
present in the [+/+] memory category in the genome- 
wide RNA-Seq data (Additional file 4: Table S4). The 
[-/-] class contained 310 genes with decreased tran- 
script levels in SI and further reduced levels in S3 



(Table 1; Figure IB). The [-/+] class contained 434 
memory genes that were down-regulated in SI, but 
produced significantly higher transcript levels in S3 
(Figure 1C); conversely, 857 [+/-] memory genes, were in- 
duced in SI but had lower levels in S3 (Figure ID). The 
latter two categories revise' their transcriptional behavior 
in a subsequent stress: after robustly responding in SI, 
these genes show weaker/no responses in S3, producing 
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Figure 1 Transcription memory response patterns of dehydration stress responding genes. Three genes from each memory category 
were randomly chosen from the RNA-Seq datasets to illustrate the diverse transcriptional response patterns. A) [+/+] memory genes; B) [-/-] memory 
genes; C) [-/+] memory genes; D) [+/-] memory genes; E) [+/=] non-memory genes; F) [-/=] non-memory genes. Transcript levels measured by 
real-time qRT-PCR are under initial pre-stressed watered (W) conditions, after the first (S1), and the third (S3) exposures to dehydration stress. 
Measurements were performed on three biological samples, each with three technical replicates. Bars are the mean ± s.e.m., n = 3 replicates. For 
interpretation of the response signs see text. 



transcripts at levels close to their initial watered (pre- 
stressed) levels. These genes are referred to as revised re- 
sponse' memory genes. The distribution of transcript 
levels for the dehydration stress-responsive memory genes 
in SI and S3 illustrating the four memory categories is 
shown (Figure 2A, B). The two non-memory gene classes 
([+/=] and [-/=]) consistently provide transcriptional re- 
sponses in S3 similar to those in SI (Figure IE, F). 

Functional distribution of Arabidopsis memory genes 

The large number of the dehydration stress-memory re- 
sponsive genes and the diversity of their transcriptional re- 
sponses raised questions of their functions and of whether 
there was a preferential association of particular cellular 
functions with any of the four memory types. To address 
these questions, genes displaying transcriptional memory 
were analyzed for biological function according to their 
Gene Ontology (GO) classification (Table 2; Additional file 5: 
Table S5). 



Genes involved in ABA/abiotic stresses responses are 
signature for the [+/+] memory genes 

Genes implicated in responses to salt, salinity, cold/heat 
acclimation, and abscisic acid (ABA) constitute about a 
quarter of the [+/+] memory genes. LEA genes represent 
3% among them (Table 2; Additional file 5: Table S5; 
Figure 3). 

Ribosome/protein synthesis, chloroplast, and thylakoid 
membrane associated functions are signature for [-/-] 
memory genes 

Memory genes encoding ribosomal, as well as chloroplast 
and photosynthesis proteins are the major constituents of 
the [-/-] memory subgroup (Additional file 4: Table S4; 
Additional file 5: Table S5; Figure 3). Encoded proteins are 
involved in ribosome structure and amino acid biosyn- 
thesis, as well as in photosynthesis, in the light harvesting 
complex (photosystem II), and in responses to light (red/far 
red light, in particular). The largest number of thylakoid 
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Figure 2 Distribution of dehydration stress responding genes in A. thaliana in SI and S3. Left panel). Transcript levels from dehydration 
stress responding genes that are up-regulated or down-regulated during SI (color key at the bottom) plotted by the log2 of their SI levels along 
the x-axis, and by log2 of their Si/watered ratio along the y-axis. Right panel) Transcript levels of dehydration stress memory genes: as in (A) except 
the y-axis is the log2 of the S3/watered ratio. The clustering of the four colors illustrates the distribution of the four distinct memory response 
types: revised response [+/-] and [-/+] memory genes clustering closer to their pre-stressed (W) levels, while the [+/+] and [-/-] increasing 
separation from these levels. 



membrane-associated proteins is encoded by [-/-] mem- 
ory genes. 

Functions encoded by [-/+] revised response memory 
genes 

There are no particular types of cellular functions specific- 
ally enriched within the [-/+] group of memory genes, al- 
though the largest percentages of genes encode chloroplast 
and membrane (plasma, organellas and thylakoid mem- 
brane) related functions. Many of these are for electron 



transport, photosystem II assembly, chloroplast organization 
and re-location, functions similar to those encoded by 
[-/-] memory genes (Table 2; Additional file 5: Table S5; 
Figure 3). 

Dehydration stress [+/-] memory genes are shared with 
multiple signaling pathways 

The signature [+/-] dehydration memory genes are regu- 
lated by multiple signaling pathways, including the ABA, 
ethylene, auxin (IAA), gibberrellic acid (GA), jasmonic 




Figure 3 Distribution of dehydration stress memory genes in A. thaliana according to GO functions. The functional GO categories 
represented by the more abundant (5% or more of the genes within a memory response type) GO groups of responding genes. The four 
memory and two non-memory types are shown in a Venn diagram format: A) the two memory and one non-memory types up-regulated in S1; 
B) the two memory and one non-memory types down-regulated in S1. 
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acid (JA), and salicylic acid (SA) pathways (Table 2; 
Additional file 5: Table S5; Figure 3). A prominent func- 
tional group represented by [+/-] memory genes encodes 
proteins associated with membranes. Although membrane 
(plasma and organellar) -related genes are highly repre- 
sented within the entire drought-responding fraction 
(including the non-memory and all four memory types), 
genes for tonoplast intrinsic proteins regulating water 
transport (GAMMA-TIP, PIP2B, TMP-C, PIP2A, TIP2, 
PIP1A, RD28) and the inward K + channel proteins 
(KAT1 and KAT2) display, exclusively, [+/-] memory 
responses (Table 2; Additional file 5: Table S5). 

Dehydration stress memory of transcription factor genes 

Among the -1500 Arabidopsis genes encoding tran- 
scription factors (TFs) [16], members of about two 
dozen families have been implicated in responses to 
drought [17-19]. Among these families, transcriptional 



memory behavior was found for members of ten of the 
families (Additional file 6: Table S6). Five families (AP2/ 
ERF, bHLH, Homeodomain-Zip, MYB/Myb-like, and 
ZF) are represented in all memory categories; GRAS and 
HSF members display only [+/+] and [+/-] memory, and 
the CAAT family is represented only by [+/+] memory 
genes. TF genes are least represented in the [-/-] mem- 
ory subgroup (Table 2; Additional file 6: Table S6). 

[+/+] memory genes for MYB2, MYB112, MYB13, 
MYB47, ATHB-7, ATHB-14, and ABF2 specifically clus- 
ter with ABA/abiotic stress responding genes but are 
rarely found among other hormonally responsive path- 
ways (Figure 4A); six putative TF genes display [-/-] 
memory type responses and only three are shared with 
other response pathways (Figure 4B). TFs with [+/-] or 
[-/+] revised response memory behavior are broadly 
shared with multiple abiotic and hormone response net- 
works, outlining distinct clusters (Figure 4C, D). 




ATMYBL2 

MYB51 

MYB28 

WRKY60 

SGS3 

AT5G54470 
SCRM2 
PIL5 
GNC 

AT2G31380 

AT4G38960 

ATHB-2 

EPR1 

DREB1A 



Figure 4 Dehydration stress memory genes encoding TFs. Hierarchical clustering of dehydration stress memory genes encoding transcription 
factors implicated in the crosstalk between various abiotic stress and hormone regulated pathways. Presence (red, value 1) or non-presence (black, 
value 0) of annotation in stress and hormone regulated pathways is determined for each transcription factor, and then clustered along both rows 
(genes) and columns (annotation) via unweighted average Euclidean distance. A) TF genes displaying [+/+] memory; B) [-/-]; C [+/-]; D [-/+]. 
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Discussion 

Biological relevance of the transcription memory 
behavior 

Plants experience more dehydration stress during the 
day when transpiration exceeds the rate of water uptake 
and fully, or partially, recover during the night when the 
rate of water uptake exceeds the greatly diminished loss 
of water from leaves [20]. This creates a diurnal cycle of 
oscillating leaf water potential that is probably most pro- 
nounced during periods of moderate drought conditions 
[21,22]. The repetitive dehydration stress system employed 
in this study has similarities to this natural diurnal stress. 
We hypothesize that dehydration stress memory helps 
plants to prepare for the next days stress if they were 
stressed the day before, despite alleviated stress signals 
during the night [9]. 

Genes displaying transcriptional memory are among 
the genes responding to the first dehydration stress. Al- 
tering their expression levels in subsequent stresses, pre- 
sumably, allows the plant to finely tune its responses to 
ongoing/recurring dehydration stress. The possible bio- 
logical relevance of the genes displaying transcriptional 
memory is considered in the context of four overlapping 
strategies generally employed by a plant to improve its 
stress tolerance and/ or survival: 1) increased synthesis of 
membrane protecting, damage -repairing, and detoxifying 
functions; 2) coordinating photosynthesis and growth 
under repetitive stress; 3) re-adjusting osmotic and ionic 
equilibrium to maintain homeostasis [23]; and 4) re- 
adjusting interactions between dehydration and other 
stress/hormone regulated pathways [24,25] . 

Increased synthesis of protective functions 

In general, the proteins encoded by [+/+] memory genes 
are involved in cell protective roles. These include heat 
shock proteins and chaperons, proteins involved in the 
repair of stress -induced damages to membranes, in 
modifying membrane lipid composition, in regulating 
membrane fluidity and permeability to toxic ions, dehy- 
drins (RD29B, RAB18, LTi30), and lipid transfer proteins 
(LTP2, LTP3, LTP4) [26,27]. LEA (late embryogenesis 
abundant) proteins functioning as molecular chaperones 
to maintain membrane structures, ion balance and 
homeostasis [28,29] and enzymes for the synthesis of 
isoleucine (toxin degradation), serine (redox responses), 
and proline (an osmolyte) are encoded by [+/+] memory 
genes. Thereby, [+/+] memory genes ensure elevated 
synthesis of factors critical for cell survival under mul- 
tiple dehydration stresses (Table 2; Figure 3). 

Abscisic acid (ABA) is the key mediator of dehydration 
signaling and is involved in responses to other abiotic 
stress response pathways [30,31]. ABA responsive genes 
that respond to dehydration are often points of overlap 
in various abiotic stress pathways [32,33] (Table 2). The 



[+/+] memory and [+/=] non-memory genes, generally, 
encode proteins with similar protective functions. How- 
ever, the memory component, presumably, allows the 
plant to optimize its stress tolerance within the context 
of protection. Only a few [-/-] memory genes and a 
smaller number of down-regulated [-/=] genes are 
shared between dehydration/ABA and other abiotic 
stress responding systems (Table 2; Additional file 5: 
Table S5; Figure 3). 

Photosynthesis and growth during repeated stresses 

More than 40% of the [-/-] memory genes encode ribo- 
some/protein synthesis, chloroplast, and thylakoid mem- 
brane related proteins indicating repetitive stress results in 
a stronger attenuation of the pathways for photosystem II 
electron transport, small RuBisCo subunit and sugar bio- 
synthesis, as well as for ribosome assembly, protein trans- 
lation, DNA replication and histones (Additional file 4: 
Table S4; Additional file 5: Table S5; Figure 3). None of 
the [+/+] memory genes encodes a thylakoid membrane 
function. Of note, 12 [+/+] memory genes associated with 
chloroplasts encode functions that are different, even op- 
posite to, the functions encoded by the [-/-] memory 
genes. Chlorophyll and toxin catabolic enzymes, enzymes 
for the synthesis of isoleucine (toxin degradation), serine 
(redox responses), and proline (an osmolyte) are encoded 
by [+/+] memory genes. Photosynthetic and starch synthe- 
sizing activities, in addition to proteins involved in chloro- 
plast organization and relocation, are encoded by [-/-] 
memory genes. Down-regulated [-/=] non-memory genes 
encode ribosome- and chloroplast-associated proteins as 
well (Table 2; Additional file 5: Table S3; Figure 3). The 
functional implications of the [-/-] transcriptional re- 
sponse patterns (decreased protein synthesis, cell growth 
and photosynthesis) and the protective functions by the 
[+/+] memory genes are in agreement with known re- 
sponses of plants enduring dehydration stress [34-37]. Of 
note, on the background of a rather even distribution of 
GO functions encoded by the late responding [=/+] and 
[=/-] genes, the fraction of genes encoding chloroplast and 
thylakoid membrane proteins constitute ~ 20% among the 
[=/-] genes (Table 2; Additional file 5: Table S5). 

Re-adjusting cellular homeostasis 

The [-/+] and [+/-] revised memory genes respond ro- 
bustly to the initial stress, but in a subsequent exposure 
provide responses at levels closer to their non-stressed 
levels. Presumably, these genes help in restoring homeo- 
stasis as the plant adjusts to dehydration stress. Thus, 
[-/+] memory genes for chloroplast organization, reloca- 
tion, photosynthetic, light harvesting, and metabolic func- 
tions are similar to the functions encoded by the [-/-] 
memory genes (Additional file 5: Table S5; Figure 3). 
However, by producing significantly more transcripts in 
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S3 than in Sl( [-/+]), or returning to pre-stressed (W) 
levels in S3 ([+/-]), these genes support the ability of 
plants to restore photosynthesis after temporary water 
withdrawal [38]. Presumably, the [-/+] genes contribute 
to an initial decrease in chloroplast/photosynthetic activity 
in SI, but by reverting to pre-stressed transcription levels 
in S3, contribute to a subsequent re-adjustment. The 
[-/+] memory genes GLK1 and GLK2, implicated in regu- 
lating the photosynthetic apparatus and chloroplast devel- 
opment in a cell-autonomous manner [39,40] are candidate 
genes for a role in such an adjustment. 

Within the same GO category, [-/+] and [+/-] mem- 
ory genes may encode opposing biochemical activities. 
For example, a number of the [-/+] memory genes en- 
code chlorophyll and carbohydrate (starch) synthesizing 
enzymes, while chlorophyll and starch degrading activ- 
ities are encoded by [+/-] memory genes (Additional file 5: 
Table S5). Juxtaposing these functions would suggest 
complementary roles during both the initial and re- 
peated stresses. Acting together, these genes are likely to 
decrease overall chlorophyll and starch biosynthesis 
during the first stress but to partially restore it during 
subsequent stresses. These complementary transcrip- 
tional patterns and functional roles support the fine co- 
ordination occurring between metabolic and energy 
adjustments during adaptation to drought [41]. 

Membrane-associated memory genes present another 
paradigm of biological relevance for the dehydration 
memory genes. As dehydration may damage membrane 
integrity, it is not surprising that a large number of 
genes (more than 750 of memory and non-memory) en- 
code functions related to membranes; more than 260 of 
the late responding genes encode membrane-associated 
proteins (Table 2; Additional file 5: Table S5; Figure 3). 
Proteins regulating osmotic pressure, water balance, and 
wall modifications have been implicated in plants' stress 
responses and environmental adaptation [42]. These 
functions are encoded by genes of all four-memory re- 
sponse types, suggesting that after providing an up/down 
transcriptional response in SI, a large number of Arabi- 
dopsis genes alter their transcript outputs (further up- 
or down-regulating or reversing towards pre-stressed 
levels) in a subsequent exposure to more finely modify 
cell membrane structure, ion balance and homeostasis. 

Dehydration stress memory genes in crosstalk with other 
response pathways 

Individual genes that are regulated by multiple different 
stimuli represent overlapping points of these plants sig- 
naling networks [43-48]. The subset of the dehydration 
stress memory genes that are at these convergent points 
revealed new aspects of plants' stress responses. For ex- 
ample the dehydration/JA pathways share the largest 
number of both memory and non-memory genes co- 



regulated by these pathways, consistent with reported 
crosstalk between jasmonic acid (JA) and ABA regulated 
signaling networks during dehydration [45,47,49], salt 
[50], and cold stresses [51]. More than 200 JA-responsive 
genes are up-regulated in SI (including memory and non- 
memory genes) and about 80 are down-regulated (Table 2; 
Additional file 5: Table S5). In S3 stress, however, 89 [+/=] 
and 46 [-/=] genes continued to provide similar transcript 
levels, seven produced higher ([+/+]), and 16 lower 
{[-/-])> transcript levels, but 140 ([+/-]) genes did not re- 
spond in the subsequent stress. Thereby, in addition to a 
finer tuning, the altered transcript levels from the memory 
genes suggest an altered crosstalk between these pathways 
in S3 compared to SI. 

Non-memory ([+/=] and [-/=]) and unidirectional 
memory ([+/+] and [-/-]) dehydration stress responding 
genes respond every time the plant experiences dehydra- 
tion stress. Concurrent with amplified responses from 
the [+/+] and [-/-] memory genes, however, engage- 
ment of the [+/-] and [-/+] genes lessens in S3, suggest- 
ing the nature and the dynamics of the interactions 
between dehydration and other signaling pathways are 
different during repeated exposures than the interactions 
occurring in a single exposure. 

Lastly, we note that TF encoded by memory genes can 
be critical for the expression of dependent genes in SI, 
but that they do not necessarily determine the memory 
behavior of regulated genes in S3. For example, the 
[+/-] memory of the MYC2 gene, identified as a master 
regulator of the crosstalk between the ABA, SA, GA, JA, 
and auxin signaling pathways [52-54], correlates with the 
[+/-] memory of a large number of MYC2-dependent 
genes common for these signaling pathways (Figure 4C; 
Additional file 6: Table S6). However, the signature gene, 
RD22, directly regulated by MYC2 [55,56] is a non- 
memory gene that is up-regulated in both SI and in S3 
(Additional file 4: Table S4; Additional file 6: Table S6) 
despite the fairly low levels of MYC2 transcripts in S3. 
Evidently, another TF activates RD22 in S3. The results 
suggest that diverse gene-specific mechanisms are in- 
volved in regulating the behavior of dehydration stress 
responding genes and that the memory behavior of in- 
dividual TFs alone does not necessarily determine or 
predict the memory or non-memory behavior of their 
targets. 

Conclusion 

The genome-wide response of Arabidopsis genes to de- 
hydration stress revealed the existence of four distinct 
transcriptional memory response patterns. By altering 
transcript levels, and presumably the levels of encoded 
proteins, memory genes are likely to alter the cellular re- 
sponses and the crosstalk between overlapping pathways. 
Adjustments of expression of memory genes, together 
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with the consistent responses from the non-memory genes, 
allow the plant to optimize its responses and the interac- 
tions between various signaling systems. Transcriptional 
memory, like defense gene priming, can provide the bene- 
fits of a more robust or modified stress response while re- 
ducing the costs of the state of preparedness [25]. The 
behavior of transcriptional memory genes adds a new di- 
mension to our understanding of plants' responses to de- 
hydration stress and to current models for interactions 
between different signaling systems. Revealing the mo- 
lecular mechanisms of transcriptional memory responses 
may be critical for understanding how plants' adapt to 
changing environments and is emerging as a new area in 
plant abiotic and biotic stress response research. 
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